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By exploring a recent model [Palmeri, J., M. Manghi, and N. Destainville. 2007. Phys. Rev. 
Lett. 99:088103] where DNA bending elasticity, described by the wormlike chain model, is coupled 
to base-pair denaturation, we demonstrate that small denaturation bubbles lead to anomalies in 
the flexibility of DNA at the nanometric scale, when confined in two dimensions (2D), as reported 
in atomic force microscopy (AFM) experiments [Wiggins, P. A., et al. 2006. Nature Nanotech. 
1:137-141]. Our model yields very good fits to experimental data and quantitative predictions 
that can be tested experimentally. Although such anomalies exist when DNA fluctuates freely in 
three dimensions (3D), they are too weak to be detected. Interactions between bases in the helical 
double-stranded DNA are modified by electrostatic adsorption on a 2D substrate, which facilitates 
local denaturation. This work reconciles the apparent discrepancy between observed 2D and 3D 
DNA elastic properties and points out that conclusions about the 3D properties of DNA (and its 
companion proteins and enzymes) do not directly follow from 2D experiments by AFM. 
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Introduction 

Whereas traditional bulk experiments provide average behaviors of dominant sub-populations, new methods exist 
that address DNA mechanical properties at the single-molecule level [T1[51I3]. Observations by AFM of double-stranded 
DNA (dsDNA) adsorbed on a 2D substrate IHH] have recently allowed a direct quantification of the distribution, p{9), 
of bending angles 6* O |7] . This led to the unexpected observation of an over-abundance of large 9 ^ , with respect to 
the Worm-Like Chain (WLC) model, at very short range (« 5 nm, much less than the persistence length « 50 nm). 
These observations suggest that, even in the absence of any bending constraints, non-linearities, such as kinks where 




3D in water , 2D in air 



FIG. 1: Sketch of a dsDNA segment solvated in water (left) with its sodium counterion cloud (the phosphate groups of the 
DNA backbone are negatively charged); and in air (right), electrostatically adsorbed on a mica substrate forming an ionic 
crystal via magnesium ion bridges between the DNA and the negatively charged substrate. The parameters associated with 
the hydrogen bonding of bps and the stacking of adjacent bases are therefore significantly modified. 

DNA is locally unstacked [3] or small denaturation bubbles, are excited solely by thermal fluctuations with a high 
enough probability to be observable at room temperature (Tr — 298.15 K). These findings cast some doubt upon the 
adequacy of the WLC model traditionally adopted in 3D [10 . In this respect, Cloutier and Widom [TT] have observed 
that short dsDNA, about 100 base-pairs (bp) long, formed looped complexes in 3D with a much higher probability 
than expected, which was attributed to partial denaturation |12j . However, these findings have been questioned by 
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new experiments that pointed out a flaw in the experimental procedure [T3] and showed that short-DNA cyclization 
data were accurately fitted by the WLC model, without invoking kinks. A recent study based on flow experiments 
draws similar conclusions [13]. These converging elements are supported by all-atom numerical simulations [HI US] 
suggesting that kinks arc not excited by thermal fluctuations with any measurable probability in unconstrained DNA 
fluctuating freely in solution. 

Apart from 2D confinement, what is the difference between both types of experiments? Figure [l] shows a sketch 
of DNA fluctuating in solution or adsorbed on a mica surface as in AFM experiments [5J IHI Ej. These experiments 
are carried out in air (the solvent is dried) and DNA is electrostatically adsorbed using magnesium ions, forming an 
"ionic crystal" with the charged substrate. DNA electrostatics are thus expected to be strongly affected as compared 
to DNA in water, hence hydrogen-bonding energies between two complementary bps and stacking energies between 
adjacent base aromatic rings are substantially modified. 

Recently, we have proposed a solvable model where bending elasticity is intrinsically coupled to bp melting |16l 117] 
in contrast to older approaches for which bending is not explicitly included [TH| [H] . Single-stranded DNA (ssDNA) 
being two orders of magnitude more flexible than dsDNA, this coupling must be taken into account because local 
denaturation strongly increases flexibility. Here, we argue that in 2D the modiflcation of the above denaturation 
parameters (bonding and stacking energies), due to adsorption, increases the probability of bp opening, which lowers 
in turn the bending rigidity. This explanation reconciles the apparent discrepancy between 3D and 2D experiments. 



Theory 
Model background 

We model dsDNA as a chain of N bps «(!<«< N) possessing two degrees of freedom [inillT] : an Ising variable, 
CTi, set to (resp. —1) when the bp is unbroken ([/) (resp. broken, B); in addition to this internal variable, an 
external one, the unit vector t^, sets the spatial orientation of the monomer. The Hamiltonian couples explicitly the 
ai and t^: 

JV-l N-l N 

H[(Ji,ti]^^ K{ai,ai+i){l-ti+i-ti)-J^ai+ia,-iJ^a,. (1) 

i=l i—1 i=l 

The bending rigidity of the joint between bps i and i + I, K{<Ti, <Ti+i), takes different values according to the internal 
state of the two neighboring bps. We denote kjj = k(1,1), kb = k(— 1,— 1) and kjjb = ^(1,— 1) = k(— 1,1). The 
Ising parameters J and /i have the following physical meanings: J is the destacking energy (energetic cost to unstack 
two consecutive aromatic rings); and 2// is the energy difference per bp between open and closed states. 

This discrete WLC model coupled to an Ising one can be completely solved using a transfer matrix approach [ini [T7] . 
Calculating the partition function amounts to solving a spinor eigenvalue problem (formally equivalent to a quantum 
rigid rotator). In 3D, the orthogonal eigenstates, denoted by l^'/^mir), are indexed by three "quantum numbers": 
? = 0, 1, . . . , oo and m = —I, . . . ,1 are the usual azimuthal and magnetic quantum numbers associated with the spatial 
orientation of and r = ± is related to the "bonding" and "anti-bonding" bp states (as for the one-dimensional Ising 
model or the -ff^ covalent bond). When projecting the eigenstates onto the real space basis lafl), with a a bp state and 
{0,(f) = fi, the two spherical angles defining t, one gets {a^\"i'i^m;T) — '4'i,m{^){<7\l,T). The ipi^mi^) = V^yi,m.{^), 
proportional to the spherical harmonics, are the eigenvectors of the pure chain model (i.e. when all k are set equal). 
The eigenvalues A^^t- are degenerate in m and can be expressed in terms of modified Bessel functions of the first kind 

{i^ = I + ^) [20] (see Ref. [17] for the expressions for the \l, r)). We have (/, t'\1, t) = Srr' , but (?', r'|/, r) ^ SwStt', 
because if I ^ I' , the matrix element is between states of different rotational symmetry. This is why our coupled 
model is not the trivial direct product of both the Ising and discrete WLC models. 

The previous exact solution can also be found when the chain is confined to 2D, as already stated by one of us in 
Ref. [21], for example when DNA is adsorbed on a substrate at thermodynamical equilibrium [7]. The spherical angles 
{0,(p) become a single polar angle 9 g (— 7r,7r]; the spherical harmonics i'l^mi&T'p) become the simpler il>n{0) = e™^, 
with n integer; the two-dimensional analogues of the eigenvalues are denoted by \n,T and the eigenvectors by |n, r) |21j . 

In the model as presented here we do not take into account additional DNA degrees of freedom, such as torsion 
or stretching. Although we have recently demonstrated that it is possible to do so in the context of thermal denatu- 
ration |22l , the additional mathematical complications of taking them into account in the calculation of the bending 
angle distribution would tend to obscure the basic physical mechanism leading to the onset of non-linear effective 
bending elasticity and is therefore not warranted here. 
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Short-distance chain statistics in 3D and 2D 

In order to compute the probability distribution p(ti ■ t^+r) of finding the polymer with a given relative orientation 
between bps i and i + r, we introduce the partial partition function, Z(zi, Zi+r), where all degrees of freedom are 
integrated out except the projections on the z axis of and t^+r, which are set to Zi and (both S [—1,1]): 

/ N r \ 

Z{z„Z^+r) = X! n / ^ 5{cOSei- Zi)5{cOSei+r - Zi+r) 

{^,=±1} Vj=i-^ j 
X {V\a^n{)\{{a,nj\P\a,+inj+{){ai,nN\V), (2) 

3 

where P is the transfer matrix and \V) the boundary vector [TB]. The complete calculation from Eq. [2] of p(s) = 
47rZ(l, s)/Z, where s = ti ■ t^+r = cosO, 9 is the bending angle between two monomers separated by a distance r, 
and Z is the full partition function, is given in the Supplementary Information, B. It uses the decomposition of P on 
the eigenbasis \^i^m;T)- We have checked that boundary effects are negligible at Tr as soon as i is larger than a few 
unities. We thus give the final result for p{s) in the limit of long DNA when the internal segment [«, z + r] is far from 
both chain ends (i.e. for iV — > 00 and i — s- 00): 

°o 97 I 1 

v{-^)^Y.-^Pi{s)Y.(Q,+%r)'e-^'^l^^ (3) 

where Pi{s) is a Legendre polynomial [20]. Eq. [3]is a sufficient approximation of Eq. SI. 12 for fitting purposes. This 
expression reveals the role of infinitely many tangent-tangent correlation lengths, ^f^ = l/ln(Ao.+/A;,r). At Tr, the 
persistence length, ~ 150 bp, coincides with the dominant correlation length [T7]. 

The same calculation holds in 2D. We find the following probability distribution (Supplementary Information, C) 

00 

n— 1 T—± 

where ^ = l/ln(Ao.+/A„ are also the tangent-tangent correlation lengths associated with 2D eigenmodes \n,T) 
with eigenvalues Xn.r- For the numerical calculation of infinite series such as Eq. |3]or Eq. ^ the sum is performed 
up to order 100 (a higher cutoff has been checked not to change numerical values). 

At room temperature, Tr, one observes below (see also Fig. [2^,c) that, for 9 smaller than a threshold 9c, p{s) and 
p(9) coincide with the discrete WLC model probability distribution, pdwlCi which is the simplified version of Eq. [3] 
or Eq. |4]when no denaturation bubbles appear (formally all n equal to ku)'- 



-2/H-l 

PDWLC(S) = 2^ -^—Pl[S) 
1=0 

PDWhciO) = E cos (n6') 



71 — — OC' 



/n(/3^) 

/o(/3'^) 



(5) 
(6) 



in 3D and 2D respectively (dotted lines in Figs. [2^, c), with (3 — {kgT)^^. In the Gaussian spin- wave approximation 
(GSW), (3k ^ 1, valid here, the discrete WLC model leads to a quadratic dependence in 9. Indeed, in this case, 
1 iPK)Y — Ii^i{f3K/r). One ends up with the probability distribution for a single joint of effective bending modulus 

k/t, andpDWLC — Pgsw — \/ (3K/{2'Kr) exp [— /3K0^/(2r)] in 2D (see Supplementary Information, E). This implies that 
the free energy required to bend the polymer by an angle 9 is quadratic, F{9,r) — K9'^/(2r). In this approximation, 
the bending rigidity k and the persistence length are related through — 2/3k in 2D and (^p = (3k in 3D [23j . 

Results 

We first examine the distribution p{s) = p(tj • t^^.^) in 3D. Whereas it is dominated at large r by the largest 
persistence length ~ 150 bp and is well described by the WLC model, this is not true at short r and large 9. 
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FIG. 2: Theoretical predictions of DNA elastic properties in two and three dimensions, (a) Logarithm of the probability 
distribution p{cos9) — p(s) in 3D (Eq. |3] full lines) for different values of r = 5, 15, 25 bp (from left to right) compared to the 
WLC model (dotted lines). One bp length is a = 0.34 nm. The Ising and elastic parameter values (in units of ksTji) come from 
fits to earlier experiments [T6j: k,u ~ k.ub ~ 147; kb ~ 5.54; fi = 1.7977; J = 3.6674. The probability distribution p{6) is given 
by p{6) = sin0p(cosS), because ds = sin 9 d6. (b) Logarithm of the probability distributions ^(6^) in 2D. Symbols represent 
experimental data taken from [7] , whereas the curves are now our best fits, from Eq. [4] The curvilinear distances between 
monomers in [7], namely 5, 10 and 30 nm, correspond respectively to r = 15, 29 and 88 bp. The value hb ~ 5.54 (in units of 
fcsTfi) comes from Ref. [TB] and ku = 160.82 comes from fitting the r = 88 bp set of data by a pure WLC model, as in [7| 
(because for such large r, the Gaussian character is restored). The remaining parameters {kub, J,tA a-re fitted. One possible 
parameter set is {kub , J, f^) =(20.97,1.3173,1.6685) (see Supplementary Information, D). Dotted line shown the predictions of 
the WLC model, for comparison, (c) Logarithm of the probability distribution p{9) in 2D. Parameter values are coming from 
fits (see Panel b), and r = 5, 15, 25 bp (from top to bottom, full lines). Dotted line shown the predictions of the WLC model 
and dashed lines show the same profiles when kb = 0. (d) Average excess chain melting AMb{6) in 2D. Same parameter 
values as in Panel b. From left to right, r — 5, 15, 25 bp. The elasticity is linear until a threshold 9c oc ^/r, where excessive 
bending induces bp melting. 



Figure [2^ displays the probability density p{s), s — ti ■ t^+r = cos 9, for realistic parameters [THIIIT]. At Tr, for 9 
smaller than a thrshold 9c, p{s) coincides with the discrete WLC model distribution, pwLc(s) (Eq. |5|, the simplified 
version of Eq. |3]when no denaturation bubbles appear. For 9 > 9^, the plot becomes non-quadratic because of partial 
DNA denaturation. The threshold 9c is estimated by equating the energetic cost of bending the polymer by an angle 
9 in its unmelted state, F{9, r) — Ku9l/{2r), with the free-energy cost of nucleating a single denaturation bubble (of 
one bp), denoted by AGb, which is AGs — 17 fcsT in 3D [17]. Using this scaling argument, we find 



2 AGs 

r (7) 
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which gives a good estimate of the observed thresholds (Figure [2p) . The anomalies (or non-linearities) appear for 
larger and larger values of 9 when r grows, and are inexistent in the plots of p{s) as soon as r > 50 bp, i.e. at 
length-scales larger than 15 nm, thus recovering standard Gaussian behavior. Indeed, setting 6c — tt in Eq.|7] yields 
the upper limit, rmax — 50 bp, as observed in the plots. This also explains why cyclization experiments with r > 50 bp 
are correctly described by the WLC model [13 . For r < 50 bp, this local melting effect is extremely weak, occuring 
with a probability Jg p{cos9) sin^^d^ « 10^ for r > 5. 

The situation is very different when DNA is confined in 2D. It has been demonstrated in experiments that DNA 
is in 2D thermodynamical equilibrium [SJ |7] . This is the reason why our statistical mechanical model applies and in 
the large N limit, the probability distribution p{6) is given by Eq. |4] Plots are provided in Figs. ^p,c for realistic 
parameter values. At large enough angles, one also sees deviations from the WLC behavior, appearing as soon as 
p{d) ~ 0.01 rad~^, a now measurable value [7]. 

We fit 2D experimental data in Fig. [2p, using Eq. [4] with kub, J and fi as fitting parameters (Supplementary 
Information, D). The fits are good over the whole 9 range. For the best-fit parameter sets, the fraction of melted 
bps for unconstrained DNA is then larger than 0.1% at Tr, two orders of magnitude higher than in 3D [16]. The 
predicted melting temperature, T^, and transition width, both on the order of 600 K, are also much higher than 
their 3D analogues. Despite the high value for in 2D, the large transition width leads, with respect to 3D, to 
non-negligible bubble nucleation, even at Tr. In other words the loop initiation factor [TH], a = e"^"^"/*-'^^^ « lO^'^ 
where Jq is the renormalized destacking parameter |17j . is increased by several orders of magnitudes with respect to 
3D [24j. The same argument as in 3D leads to r^^yr ~ 120 bp in 2D, after modifying AGb — 6.6 kgT according to 
our fitted parameters. Furthermore, we display in Fig. |2jl the average excess of melted bps when • t,;+r = cos 6* is 
fixed, as compared to an unconstrained DNA (see Appendix). As anticipated, the deviation from the WLC behavior 
at 9c coincides with the appearance of melted bps making the polymer more flexible. 

Discussion 

How can the apparent discrepancy between 2D and 3D parameter values be explained? Not by the fact that the 
DNA used in 2D experiments are heteropolymers, whereas the values derived in 3D come from poly(dA)-poly(dT) 
homopolymers [11]. Indeed, even for the most robust poly(dG)-poly(dC), Tm — 360 K in solution. A simple and 
straightforward explanation for the discrepancy in parameter values is related to the change in the DNA electrostatic 
energy when it is solvated in water (3D) or adsorbed through magnesium (Mg^+) bridges on the mica in a dry 
environment. Indeed, it is known that slightly modifying electrostatic interactions (such as by varying the salt 
concentration) changes dramatically the denaturation profile of DNA in solution (see e.g. |25j). The energy required 
to break a bp, 2/i, and the energy to destack consecutive bps, 2 J, should also be sensitive to the change in the direct 
adsorption energy between mica and ds or ssDNA. Strong support for this mechanism comes from the experimental 
results of Wiggins et al. themselves [7 . In their Fig. S3, they present the angle distribution and cnd-to-cnd distance 
statistics for DNA adsorbed on a different-quality mica. Even though the data match to a good approximation those 
of their Fig. 3, a detailed analysis of the plots for r = 5 and 7.5 nm leads to the conclusion that the two data sets do 
not coincide, even taking into account error bars. This is an experimental indication that the substrate on which DNA 
molecules are adsorbed does indeed influence its microscopic parameters. Recent AFM experiments also testifled to 
a DNA structural modiflcation after adsorption on mica and drying |26j : poly(dG)-poly(dC) proves to shorten its 
contour length, supposedly by taking an A-DNA conformation, in contrast to poly(dA)-poly(dT) or plasmid DNA, 
both of which keep their B-DNA conformations. 

As a result, inferring the parameters fi and J from their 3D analogues is a challenging task. At the present time, 
the best strategy is certainly to flt them to experimental data. The above results are confirmed by recent accurate 
all-atom molecular dynamics simulations: Mazur has investigated in detail the short-distance angle distribution of 
3D DNA and did not find any evidence for the strong deviations from a WLC distribution found experimentally in 
2D [H]. 

Now we discuss in greater detail the role of bubble fiexibility, k b , and of cooperativity, J, by comparing our model 
to earlier ones. In the kinkable WLC model [27], kinks of vanishing rigidity can be activated by thermal fluctuations. 
This model and ours become physically equivalent in the — *■ limit: a 2 bp denaturation bubble plays the role 
of a "kink", in the sense of a thermally activated local defect without rigidity. Our microscopic vision of a kink 
thus differs from Lankas et a/.'s local unstacking one [9], but yields the same short-range mechanical properties. 
When Kb = 0, the interesting behavior of p{9) in the denatured region is destroyed: p{9) becomes flat (Fig. [2j:), 
as in |27j . and is practically insensitive to r once a kink is nucleated, because a chain segment including a kink has 
vanishing rigidity. This is the reason why Wiggins et al. appeal to a different Linear Sub-Elastic Chain (LSEC) model, 
with a phenomenological bending energy -Blsec = ^I^L which enables them to satisfactorily fit their experimental 
data |7, ,28, . In contrast to this LSEC model, our approach proposes a microscopic explanation associated with bubble 
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nucleation for the sub-harmonic behavior of p{0). Due to excess bubble formation, our model predicts deviations from 
WLC (or Gaussian) behavior as soon as r < rmax with Tmax = '^'^ ku / {2AG b) (from Eq. [7|. This expression differs 
from the LSEC model one, for which r^ax ~ Pkij. 

Setting J = with kb finite also affects the profiles by softening the transition and increasing significantly the large 
angle probabilities, by a factor greater than 10 (data not shown), which confirms the importance of cooperativity 
(when in addition kjjb — 0, the model proposed in Ref. [T2] in the context of cyclization is recovered). Neglecting J 
or Kjjb would require the use of unphysically large kb values when fitting experimental data, while worsening the fit 
quality. 

Our model is restricted to homopolymer DNA. However, a more accurate treatment should incorporate sequence 
effects by using bp dependent model parameters [29]. Considering that the heteropolymer case is difficult to treat 
theoretically, and experiments provide only an average description of bending angle probability distribution, we limit 
ourselves here to describing the anomalous behaviour using an averaged approach. If more detailed experimental 
results become available, it would be worthwhile to extend our model to treat the heterogeneous case. 

Currently, many AFM experiments explore DNA conformations and complexation between nucleic acids and pro- 
teins (see reviews [UlSniEI])- When AFM imaging is carried out on DNA [B', [71 [SH 133] or DNA/histone complexes [M] 
in order to access their statistical and dynamical properties, effects of surface interactions on DNA structure are likely 
to modify sensibly these properties. More generally, our work suggests that studying DNA/companion proteins 
interactions by AFM [ 351 ESI 123 EH] does not provide any quantitative clue to 3D complexation. 

In the cell, packaging involves wrapping DNA around positively charged histones [39 . It has been shown that 
this adsorption is mainly driven by electrostatics [ID]. Our results suggest that in this case, DNA adsorption on a 
curved charged surface (such as the histone) is likely to modify profoundly local elastic and denaturation properties 
of dsDNA. Enhanced flexibility due to denaturation is then likely to facilitate wrapping. This mechanism might also 
be important for improving the accessibility of enzymes to the single strands in local bubbles [411 142j when DNA is 
wrapped. 
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FIG. 3: Logarithm of the probability distributions p{9) in 2D (r = 15 bp), for both the LSEC model (with A = 27.2 pN.nm [7], 
dashed lines) and our theory (full lines), for increasing temperature T. From top to bottom, T = 298.15 K, 330 K, and 360 K. 
Our model predicts that increasing T enhances flexibility in a pronounced manner (note the logarithmic scale) thanks to the 
opening of bps. At 6 = 1.5 rad, an experimentally accessible value 0, the decrease of — In p{9) with increasing T is twice that 
found with the LSEC model. 



One way of validating the present model at the experimental level would be to quantify the effects of temperature, 
which can be predicted for both our coupled model and the LSEC one [55] (Fig. [3] see Supplementary Information, 
E, for LSEC formula). Our model predicts that increasing temperature enhances flexibility in a more pronounced 
manner, thanks to the opening of bps. We believe that such a deviation between the predictions of both models 
would be a credible experimental test of their respective validities. Additional tests of the quantitative difference 
between DNA properties in 3D and 2D would be to compare cyclization rates by AFM in both situations for the same 
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dsDNA strands, or to check that denaturation remains weak in 2D when approaching the 3D melting temperature, 
as predicted by our results. 
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Appendix: Bending-induced melting in 2D 

Following a calculation as in Ref. P7j, we derive the excess chain melting AM^ as a function of 9. It measures 
the average excess of melted bps in the bended chain as compared to the free, unconstrained one and is given by 
/S.Mb{6) = — ^^g^lnj>(0) (see Supplementary Information, F). The comparison of Figs. 2c and d confirms that the 
deviation from the WLC model corresponds to the appearance of melted bps that make the polymer more flexible 
at short range. An interesting feature of these calculations is the saturation of AiV/s at a finite value, even when 
r < Tinax increases. In Fig. |2}i, this value is close to 3, which means that the total excess number of denatured bps 
does not exceed 3 on average. In other words, even if r bps, or more, can in principle be melted to relax the constraint 
ti • ti_|_r = cos{9) , only a few of them actually do, since it costs more energy to melt more bases, whereas, owing to the 
small value of kb, a small denaturation bubble suffices to give the whole molecule a very small resistance to torque. 
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A. Effective Ising Hamiltonian 

Starting from the Hamiltonian H[ai,ti\, an effective Ising Hamiltonian, a function of the only, is obtained by 
integrating out the rotational degrees of freedom in the partition function, leading to a renormalized Hamiltonian, 

JV 

-^f^CTj, (1) 
1 = 1 

with renormalized parameters: Jq = J — kBT[Go{(3Ku) + Gq{[3kb) — 2Go(/3k;7b)]/4 and = —kBT[Go{(inu) — 
Go(/3kb)]/2 where (3 = {kBT)~^ and Gi){x) = x — ln(sinha;/a;) is related to the bending free energy of a single joint; 
/i is not renormalized (see Ref. [17' for further details). 



E 



B. Calculations of the probability distribution p(cos 0) in 3D 

In Eq. (3) of the body of the text, we define the partial partition function where all degrees of freedom are 
integrated out, except the projections on the z axis of t,j and t^j^j., set respectively to Zj and Zi^^.. Imposing the value 
of ti ■ ti^r = s £ [—1, 1] amounts to fixing Zi = 1 (i.e. t^ — z), and Zi+r = s (see Fig. SI. [T]) and to multiplying by the 
solid angle An to restore the rotational invariance of the whole problem with respect to t^, which can actually take 
any orientation. Thus we have p(ti ■ tj+r = s) = 47rZ(l, s)/Z where 

Z={V\P^-^\V) (2) 
is the total partition function (see Methods and Ref. [17]). 




Fig. SI. 1: Examples of chain configurations counted by the same partial partition function Z{l,s) (Eq. (3) of the body of the 
text), in the case where the singled-out tangent vectors lie on the polymer ends, = ti = z, t^+r = tjv and s = z • tjv = cos 6. 

We denote by the projector on the z axis defined by H^lcrri) = S{cos9 — z)\<t^). It follows from Eq. (3) of the 
body of the text that 

Z{z,, z,+r) - {V\P'-'KP'-ii,^^^P^—'\V). (3) 
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In order to compute this quantity, we need to express in the basis where P is diagonal, namely the l^i^mir)- We 
use 

{an\i'i,rn,r)^^l,mm{cT\l,T), (4) 

where the spherical harmonics ■0/,m are defined by the Associated Legendre polynomials P/ ,„ PU| as follows: 

^Pi,mm = ^^,™(0,(^) = fc^,„P^,|,„|(cos^^)e™'^, (5) 



Below, we shall only need the Legendre polynomials Pi{x) — P/,m^o(^) [20 . Now we can compute 

-^{^l',ni',T'\flz\crn){an\i'i^^,r) (7) 
— {^l,,,n';r'Wn)S{cOSe- Z){an\^l,,n;r) (8) 



(9) 



= ]^h,mkl',rnPL\m\{z)Pi>^\rn\{z){l' ,T'\l,T)5„,,m'- (10) 

When m — m' = 0, the previous equality specializes to 

(^',,,0;r'|n.|^'/,0;r) = ^V2r+lV2FTlPi{z) Pi, (z) {I' , t' \l , t) . (11) 

If the boundary vector \V) has rotational symmetry (/ = 0), then {^i^mirlV) — ^io^mo{Oy'''\V)- It follows that 
47rZ{z,,z,+r) -EE ^-(^|0,T)A^7;Po(^0(0,r|?,rOP/(zOAL.Pz(^.+,0(^,T'|0,T'OPo(^.+.)A^7,r'(0,r''|y). 

{t.t'.t"} (=0 

(12) 

This partial partition function must be compared to the complete one 

Z = {V\ P^-i \V) = J2\{y\0,r)\'\^-' (13) 

r 

in order to get 



S^o(2/ + l)(^|0, r)A^;^(0^^ 

2E.I(^|0,r)PA„^7 



p(t,; • t,+, = S) = ^ ^|2^^^-l ^ 



because Po(l) = ^0(2;) = P/(l) = 1- One can check that the distribution is correctly normalized, J^^pis) ds — 1, 
because J^^ Pi{s) ds — 2Sifi. 

In the limit of a long DNA where the internal segment + r] is far from both chain ends, in other words when 
N —> 00 and then i 00, the previous relation becomes 



00 



p{U-U^r = s) = ^?i±lp,(.) ^ |(0,+|?,r')n^) (15) 
= 2+£'^^'(^) S |(0,+|/,r')pe-/«-. (16) 

1=1 r'=± 

This expression reveals the role of infinitely many correlation lengths, the — l/ln(Ao,+/A;_7-)- At Tr, the persis- 
tence length ~ 150 bp coincides with [TSJ[T7]. We have checked that boundary effects are indeed negligible at 
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room temperature (Tr = 298.15 K) as soon as i is larger than a few unities. Thus Eq. (16 1 is a sufficient approxi- 
mation of Eq. (14) for fitting purposes and is used in the body of the text. Once this distribution p{cos6) is known, 
the probability distribution of 9, denoted by p{0), is simply given by p{9) — p{cos 6) sin 9d9 , because s — cos9 and 
\ds\ = sin 6* \d9\. 

As a corollary, the mean value of the correlator (t^ • t^+r) can be computed in this limit: 



s p{s) ds 



f2^^ ( f\ p^s) ds) y: m+\iy) 

1=0 ^ ^-^-1 ^ r'=± 



r' = ± 



Ai,t' 



(17) 
(18) 

(19) 



/■I 2 

because / s Pi{s) ds — -Si^i. We therefore recover the result of Rcf. [17 , Eq. (100), where it is pointed out that 



only two correlation lengths remain in this correlator. 



C. Calculations of p{9) in 2D 

Using the 2D algebraic background presented in the Methods and the fact that {a9\^n;T) — e"^^ {a\n, r), the matrix 
elements of the projectors II^ in the eigenbasis become 

{^n';r'\tiz\^n;r) = l ^ (*«'.r' |'Te>^(^ " ^o) (a^|*n;r) (20) 

= ^E^^^"""'^'''("''^V)(^i'^'^) (21) 

CT = ± 

= -^e'("-"')^«(n',r'|n,T). (22) 
27r 

In addition, (^'ri:T|^) = <^rio(0, t\V) by the rotational symmetry of the boundary conditions, and comparing the partial 
partition function, now denoted by Z{9), with the full one, Z, leads to the 2D counterpart of Eq. (16 1: 

p{9r^9) ^ ^ (23) 



I- E oos{n9) E l(0,+|n,r')P (^Y (24) 



2^ n£"oo V ^0.+ 

1 + 1 g cos(n0) E (0. +1". ^? e-'-Z^S,. , (25) 
27r TT ^-^ 

n— 1 T— ± 



valid in the limit of long DNA strands, with 9^ E (— tt, tt] defined by cos 9r = t,; • ti^r and ^ = 1/ ln(Ao.+/A„_7-). 



D. Fitting p{9) in 2D 

Our first approach to fitting the 2D experimental data of Fig. 2B of the body of the text, coming from Ref. [7], 
consisted in directly using the 3D parameter values, in particular those of J and /i, in the 2D Hamiltonian. The 
so-obtained angle probabilities are very far from the experimental ones, which led us to fit the model to experiment. 



We display in Fig. 2B of the body of the text our best model fits, using Eq. (25) with kjjb, J and /i, as fitting 
parameters. These least-square fits are good over the whole 9 range. The bp length is assumed to remain a — 0.34 nm. 
Thus the curvilinear distances between monomers in [7], namely 5, 10 and 30 nm, correspond respectively to r = 15, 
29 and 88 bp. The value (in units of ksTfi) of Kg = 5.54 comes from Ref. [TB] and ku = 160.82 comes from fitting 
the r = 88 bp set of data by a pure WLC model, as in [7] (assuming that for such a large r, the Gaussian character 
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is restored thanks to the central-limit theorem; we have checked that it is indeed the case for the parameters that 
we discuss now). The three remaining parameters {KuB^J^^^) are fitted by a simulated annealing algorithm of our 
own. The best fit values appear to be highly degenerate, in the sense that a whole subset of parameters in the 
three-dimensional {kub, J, fJ-) space yields essentially the same mean-square deviation. 

This degeneracy can be related to its 3D equivalent. In Rcf. [17 , we have simplified the discussion by setting 
Hub = /"^u, the value for dsDNA, because varying kub amounts simply to changing the bare value of J in order to 
keep the melting temperature, T^, and the transition width unchanged: Jmod = J + ^^[Go(/3k;7) — Gq{Pkub)] |17j . 
However, when it comes to the angle distribution p, the value of kub might play a more fundamental role because the 
different persistence lengths depend on it via the eigenvalues A in Eq. (25 1. In 3D, we have checked that changing the 
value of KjjBi while suitably modifying J, has, in practice, little influence on p, Nr, and AMb (see definitions below), 
even at very short scales (r = 5 bp). In 2D, suitable values of kjjb range between 10 and 50 kBTn. Even if kub = 
or Kub = f^B is held fixed, the minimization with respect to J and /i alone does not give a significantly poorer fit. 
Examples of parameter sets provided by simulated annealing are, in units of fc^Tfj, (kc/b, J, m) =(20.97,1.3173,1.6685) 
or (45.10,0.8637,1.7885). Increasing the temperature by 20% does not lift the degeneracy. With these values, the 
fraction of melted bps for an unconstrained DNA varies between (ps ~ 0.1% and 0.4% at Tr. 



E. Probability distribution p{9) calculated from diverse Hamiltonians in 2D discussed in the paper 



We consider a pure elastic chain without bp melting, described by a 2D single-joint Hamiltonian H{0). Without 
any loss of generality, the probability distribution p{9) in 2D can be written as 



p{9) 



riLi /-^ de, exp [-l3H{e,)] 5{6 ~ g,,,+.) 
m=i/-. '^^»exp[-/3ff(&0] 



(26) 



where 9 = 0^,^+^ — X]}=i ^ bending angle between t^ and tj+j.- By introducing the Fourier transform of the 

5 distribution, we get 



p{9) 



1 

2^ 



due'' 



z{uj) 



where the characteristic function of the single-joint Hamiltonian is defined by 



z{uj) 



[ d9 exp[-iuj9 ~ I3H{9)]. 

J — TT 



(27) 



(28) 



We consider four cases: 



1. For a Gaussian Hamiltonian (Gaussian Spin Wave approximation, GSW), H{9) = k9'^ /2 with the approximation 
9 G (— cxD,oo), we have z{uj)/z{0) = exp[— w^/(2/3k)] and we easily get the Gaussian probability distribution in 
2D 



PGSW(^') 



f3K 

2-Kr 



exp 



2r 



(29) 



2. For the discrete wormlike chain where H{9) = k(1 — cos 9) is periodic in 6* G (— tt, tt], Eq.(27l reduces to a 
decomposition in Fourier scries of modes n with z{ijj) becoming e^^'^IniPn) and we find 



^ oo 

PDWLc(6') ^ cos(n6') 



n— — oo 



(30) 



This is equivalent to Eq. (p5| in the simple case where all the n are set equal. 



3. For the LSEC model (T) |2S], where H{9) = h.\9\ with the approximation 9 E (— oo,cx)), we have z{uj)/z{0) 



1/[1 -I- (w//?A)^] and the probability distribution 



PLSEC(^') 



2k-rpA {f3A\9\r-'^K^_^{f3A\9\) 
r(r) 



(31) 



where F is the Euler function and Ki,{x) the modified Hankel function [20 . Equation (31 1 is plotted for different 
temperature values in Fig. 3 of the body of the text. Note that in the case where 9 G (^tt, tt], the above result 
has a negligible correction to p on the order of exp(— /JAtt) ~ 5 x 10"""^° rad~^ if /3A ~ 6.8 
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4. Both the GSW and LSEC models are special cases of a generalized model where H{9) = a\6\^. With the 
approximation 6 € (—00,00), we recover the GSW model when rj = 2, a = k/2 and the LSEC model when 
rj = 1, a = A. Although the general case is harder to handle for arbitrary r, two limits, r = 1 and large r, are 
easily studied. When r = 1, 

For large r, [z(a;)/z(0)]' becomes a sharply peaked Gaussian function centered at w = 0: [z{uj) / z{0)]'' ~ 
exp[— rw^(6'^)/2]. In this limit p{9) becomes effectively Gaussian, 



P(^) - exp [ ] ^ (large r), (33) 



2nr V 2r 



1 



can 



where Keff(a,?7;/3) = l/(/3(6'2)) = (/3a)2/'T(l/77)/[/3r(3/?7)]. The effective bending rigidity, k^s = 
be written in terms of the asymptotic persistence length, the segment length, £, and the temperature. It is 
only for rj = 2 that p{0) is Gaussian at all length scales r. For strongly subharmonic models (77 well below 2), 
the large angle bending distribution can be much higher than the Gaussian model prediction for r <C ^ (see 
Fig. 3A of Wiggins et al. [7] for the 77 = 1 case). For non-Gaussian models, p{9) interpolates smoothly between 
an "anomalous" behavior [Eq. ( [32| )] and an effective Gaussian one [Eq. ( [33| ] as r increases from one past the 
persistence length, ^. For large r, even at large 9, p{9) is effectively Gaussian because it is determined mainly 
by those high entropy configurations for which almost all joint angles, 9j, are small (therefore in this case only 
the small 6 behavior of H{9) is important). 

From fitting the large r data to the Gaussian model, Wiggins et al. found (3kcs — ^jl = 54nm/2.5nm ~ 22, 
which implies that /3A = /3aLSEC = (2/3Keff)^^^ = 6.6, close to the value of 6.8 found by fitting their LSEC 
model to experiment using Monte Carlo simulations. 



F. Bending-induced melting in 2D 

Following Wiggins et al. [27j , we now derive the average bending moment Nr as a function of the deflection angle 
9. In the 2D case, 

Nr^-^lnZi9) (34) 

(in units of ksT) measures the torque perpendicular to the substrate that must be applied to two interior monomers 
separated by r bp in order to impose a deflection angle 9 between them. Examples of vs 9 plots for our model are 
given in Fig. SI. [2] Four regimes appear in these plots: 

(i) a linear one at low deflection angle 9 < 9c\ 

ill) a non-linear one for intermediate angles; 

(iii) a saturation plateau for large deflection angles and 

(iv) a decreasing one near 9 — 

In regime (i), the DNA response is linear, with a GSW bending moment = f3dF{9,r) /d9 — {£,^/r)9, determined 
only by ~ finui since melted bps are essentially inexistent. In the intermediate region (ii), a non-linear behavior 
occurs. Indeed, it becomes more favorable for the system to break bps in order to make them more flexible and 
thereby relax the high bending constraint. The plateau appearing in region (iii) shows non-zero response due to the 
finite value oi kb, contrary to the kink model [27j. As already mentioned in this Ref. when 9 approaches tt 
[regime (iv)], A^,. vanishes: the symmetry of the system through the axis defined by the vectors tj and ti^j. imposes 
that = TT is an (unstable) equilibrium point. Indeed, 9 S (— 7r,7r] in our calculations, thus larger angles are brought 
back in this interval modulo 2tt. For a given 9, the summation in Z is in fact a sum over all the 9 + 2fc7r, fc e Z, as 
illustrated in Fig.SI. [l] Since the contributions oi 9 = n and 9 = — tt cancel, A^^ vanishes at 9 ~ tt. More importantly, 
we cannot give quantitative answers for cyclization experiments when \9\ ~ 27r. 
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Fig. SI. 2: Torque Nr (top) and excess chain melting AMb (bottom) as a function of 9, in 2D. Same parameter values as in 
Fig. 2B of the body of the text. From left to right, r — 5, 15, 25 bp. The elasticity is linear until a threshold 6c oc y^, where 
excessive bending induces bp melting. 

To render these physical mechanisms more explicit, we have also displayed in Fig. SI. [2] the excess chain melting 
AMb as a function of 6 (Fig. 2D in the body of the text). It measures the average excess of melted bps in the 
bended chain as compared to the free, unconstrained one and is given by AMb = l'^P(^)- We thus confirm 

that the collapse of Nr corresponds to the proliferation of melted bps. The typical angle 9c at which bending- induced 
melting occurs is again estimated by equating the energetic cost of bending the polymer in its unmelted state, of order 
PKu0c/2r, with the free-energy cost of nucleating a single denaturation bubble (of one bp), AGs = 'iJo+2KQ + 2n |17j . 
Again, this argument that leads to 9c ~ y/r gives a good estimate of the observed threshold. It also gives the upper 
limit of r for which these non-linearities are apparent, Tmax ~ ^ku/AGb — 120 bp in 2D. 

An interesting feature of these calculations is the saturation of AMb at a finite value [regions (iii) and (iv)], even 
when r < r^ax increases. In Fig. SI. |2] this value is close to 3 and the total excess number of denatured bps does not 
exceed 3 on average. This is corroborated by the fact that the large 9 torque is independent of r, mainly due to the 
few melted bps. In other words, even if r bp, or more, can in principle be melted to relax the bending stress, only a 
few of them actually do, since it costs more energy to melt more bases, whereas, owing to the small value of kb, a 
small denaturation bubble suffices to give the whole molecule a very small resistance to torque. 

Note that iV^ is essentially unchanged between the different fitted parameter sets discussed above, and AAIb varies 
by at most 20%, in the melted region only. 



G. Bending-induced melting in 3D 

The previous calculations can be extended in a straightforward way to the 3D case, with very similar qualitative 
conclusions. The physical meaning of the bending moment is less direct because of the axial symmetry (m = 0) 
imposed in the calculation of p(cos6'). By contrast, the excess melting, AMb, is meaningful for circular DNAs, where 
9 is close to 27r, with the reserves given above. The above argument now leads to rmax ~ 50 bp. In Fig. SI. |3]one 
actually sees that for r < 50 bp, AMb saturates at 10 bp near 9 = tt. In the topical case where r is comparable 
to the chain persistence length (~ 150 bp), bending-induced melting of constrained DNA is expected to be virtually 
inexistent in 3D. As for large looped complexes, such as in Ref. [3] where r ^ 1000 bp, melting is not expected to 
stabilize or to facilitate looping either. 

In contrast, non-linear behavior can play a significant role by making DNA much more flexible when the polymer 
is highly bent, such as in short circular DNA (see [9]) or in protein-DNA complexes. Our calculations on bending- 
induced melting, AMb, give a good indication of the degree of melting following a sharp bending constraint. Note, 
however, that a complete treatment of cyclization would require imposing not only the angle 9 but also the physical 
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Fig. SI. 3: Torque Nr (top) and excess chain melting AMb (bottom) as a function of 9, in 3D. Same parameter values as in 
Fig. 2A of the body of the text. From left to right, r = 5, 15,25. The elasticity is linear until a threshold 6c oc ^/r, where 
excessive bending induces bp melting. 

distance between bps i and i + r 12\. Although the full calculation is outside the scope of the present work, imposing 
6 does already contain some of the important physical features of cyclization. It is reasonable to expect that even for 
~ TT, our calculation reproduces the correct order of magnitude for cycled DNA. For instance, our 3D predictions 
concerning excess melting, e.g. AMb ~ 10 bp, could be checked by doing UV absorbance measurements on short 
circular DNAs. 





